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RESUME. Dans le but de reduire la masse des panneaux solaires d'un satellite, la flexibility de ces 
panneaux est devenue non negligeable et ils peuvent se taper I'un contre I'autre ce qui peut provo- 
quer un endomagement de la structure. Pour prevenir cela, des cales sont fixees dans des positions 
bien choisies de la structure; ils jouent le role d'un ressort elastique unilateral. Comme consequence 
negative, la dynamique de ces panneaux devient nonlineaire. Une approximation par elements finis 
est utilisee pour resoudre les equations aux derivees partielles qui gouvernent la dynamique de la 
structure; des balayages en frequence ont ete faits numeriquement pour etudier le comportement 
dynamique. Les modes normaux non linaires sont a I'etude. 

ABSTRACT. In order to save mass of satellite solar arrays, the flexibility of the panels becomes not 
negligible and they may strike each other; this may damage the structure. To prevent this, rubber 
snubbers are mounted at well chosen points of the structure and they act as one sided linear spring; 
as a negative consequence, the dynamic of these panels becomes nonlinear. The finite element ap- 
proximation is used to solve partial differential equations governing the structural dynamic. Frequency 
sweep has been performed numerically to study the dynamic behavior. Non linear normal modes are 
under study 

MOTS-CLES : Dynamique non lineaire, Systemes differentiels, Moderation 
KEYWORDS : Nonlinear dynamics, Differential systems, Modeling 
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1. Introduction 

In articles Q], [3], 0] and |5], the dynamic of a beam system with a nonlinear contact 
force, under a periodic excitation given as an imposed acceleration form is studied both 
numerically and experimentally. When sweeping frequencies in an interval which contains 
eigen frequencies of the beam, resonance phenomena appear as well as new frequencies 
caused by the unilateral contact. 

The study of the total dynamic behavior of solar arrays in a folded position with snub- 
bers are so complicated, that to simplify, a solar array is modeled by a clamped-free 
Bernoulli beam with one-sided linear spring. This system is fixed on a shaker which has 
a vibratory motion d(t), it is given as a periodic imposed acceleration : d(t) = asmujt 
(see figureQ}. 

The present study is to simulate the behavior of a beam which strickes a snuber. In 
this first step, existing algorithms are used in order to prepare a sequence of experiments. 
The focus is on the comparison between the linear and the non-linear case. 
An original approach of an extension of normal modes to this non-linear case is under 
study ; this approach is an alternative to the one of J2] and 
The beam motion with a snubber can be modeled as : 

pSu(x,t) + EIu [lv) {x,t) = [1] 

with the boundary conditions : 

u(0,i) = d(t), d x u(0,t) = 0, 

EIuW(L,t) = Oand EIu^{L,t) = k r (d(t) - u{L.t))+. 

u(x t)+ = l l f u>0 m 

[ ■ >+ \ ifu<0 L J 

The classical Hermite cubic finite element approximation is used, it yields an ordinary 
differential system in the form : 

Mq + Kq = k r (d{t)-q n ) + e n [3] 

Where M and K are respectively the mass and stiffness assembled matrices, q is the vec- 
tor of degrees of freedom of the beam, = (itj, d x Ui), i = 1, n, where n is the size 
of M. 

The system is integrated numerically using the Scilab routine 'ODE' for stiff problems, 
package ODEPACK is called and it uses the BDF method (§]. 

The frequency sweep is performed by running the integration of ordinary differential 
systems at a chosen initial value of the frequency and saving the maximum of the displa- 
cement of all nodes in the whole time interval. The frequency is then incremented till the 
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Frequency sweep for a beam system with 3 



end of the frequency interval and the same procedure is used again. 
The maximum of displacement is plotted in the frequency domain, the simulation shows 
many resonance phenomena when approaching particular frequencies ; these particular 
frequencies are the non-linear frequencies of the system corresponding to the eigen fre- 
quencies of the linearized system (bilateral spring). Other superharmonic and subharmo- 
nic frequencies appear due to the nonlinearity of the contact. The frequency responses are 
compared with the FFT of the temporal signal coming from the direct integration of the 
differential systems. On the other hand, the calculation of the non-linear normal modes 
(NNM) [3 1 and [2] will be a theoretical way to validate the numerical simulations. 




Figure 1 — At left : solar arrays from folded to final position, at right : beam system with unilateral 
spring on a shaker to model an unfolded array 



2. Numerical results 

The Physical properties and the dimensions of the beam with a spring stiffness of 
5 x 10 5 N/m gives the following eigen-frequencies (bilateral spring) : 

- first frequency — 196.35694.ffz second frequency = 472.07584.ffz 

- third frequency — 961.52462.ffz 

The figure (j2| show the eigen-frequencies of the linearised system (bilateral spring) cal- 
culated analytically and the corresponding frequencies of the non-linear cases (unilateral 
spring) ; many other frequencies appear in the non-linear cases due to the nonlinearity 
of the contact. The superharmonics and the subharmonics frequencies are shown too, the 
superharmonics are the double and the triple (etc..) of the corresponding eigen frequency, 
the subharmonics are the half and the third (etc..) of the eigen frequencies. Many other 
combinations of these frequencies may exist but it is quit difficult to distinguish their lo- 
cation. 

A damping term in the spring will help to distinguish each mode with a particular initial 
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condition or excitation, the full article will show this technique. 

The numerical results are to be compared with experiments in preparation. 

On the other hand, the calculation of the non-linear normal modes are under development 

in order to address more precisely the frequency response of the structure and to validate 

numerical results. 

Comparison with results obtained by asymptotic expansions of the type in 6 J is in project 



Frequency sweep, unilateral and bilateral 
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Figure 2 - Frequency sweep, bilateral and unilateral spring, no pres stress or backlash, 

beam height = 10mm, kr — 500.000./V/m, a = 50m/s 2 , tf = OAs 
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